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lON-BY-ION COOLING EFFICIENCIES 

Orly Gnat^'^ and Gary J. Ferland^ 

ABSTRACT 

We present ion-by-ion cooling efficiencies for low-density gas. We use Cloudy (ver. 08.00) 
to estimate the cooling efficiencies for each ion of the first 30 elements (H— Zn) individually. 
We present results for gas temperatures between 10* and 10^ K, assuming low densities and 
optically thin conditions. When nonequilibrium ionization plays a significant role the ion- 
ization states deviate from those that obtain in collisional ionization equilibrium (CIE), and 
the local cooling efficiency at any given temperature depends on specific non-equilibrium 
ion fractions. The results presented here allow for an efficient estimate of the total cooling 
efficiency for any ionic composition. We also list the elemental cooling efficiencies assuming 
CIE conditions. These can be used to construct CIE cooling efficiencies for non-solar abun- 
dance ratios, or to estimate the cooling due to elements not included in any nonequilibrium 
computation. All the computational results are listed in convenient online tables. 

Subject headings: ISM:general - atomic processes - plasmas 



1. INTRODUCTION 

The radiative cooling efficiencies of hot (lO** — 
10® K) low density gas are important quantities in 
the study of the diffuse interstellar and intergalac- 
tic medium. They determine the thermal and dy- 
namical properties and evolution in a variety of as- 
trophysical systems, ranging from local interstellar 
clouds to shocks in intergalactic filaments. 

Computations of hot gas cooling efficiencies as- 
suming collisional ionization equilibrium (CIE) have 
been studied extensively (House 1964; Tucker 
& Gould 1966; Allen & Dupree 1969; Cox 
& Tucker 1969; Jordan 1969; Raymond et 
al. 1976; Shull & van Steenberg 1982; Gaetz 
& Salpeter 1983; Arnaud & Rothenflug 1985; 
Bochringer & Hensler 1989; Sutherland & Do- 
pita 1993; Landi & Landini 1999; Benjamin et 
al. 2001). These CIE cooling efficiencies depend 
only on the gas temperature and metallicity. How- 
ever, there are many cases in which CIE condi- 
tions do not apply. For example, nonequilibrium 
ionization is bound to occur when an initially hot 
gas cools radiatively below ~ 10^ K (Kafatos 1973; 
Shapiro & Moore 1976; Edgar & Chevalier 1986; 
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Schmutzler & Tscharnuter 1993; Sutherland & Do- 
pita 1993, and Smith et al. 1996; Gnat & Stern- 
berg 2007). Below this temperature, cooling be- 
comes rapid compared to electron-ion recombina- 
tions, and the gas at any temperature tends to 
remain "overionized" compared to gas in CIE. In 
conduction fronts surrounding evaporating clouds 
(e.g., Borkowski et al. 1990; Gnat et al. 2010), non- 
equilibrium ionization occurs when the ionization 
time is long compared to the rate of temperature 
change. In this case, the gas tends to remain unde- 
rionized compared to CIE. Non-equilibrium ioniza- 
tion also plays a role in fast radiative shock waves 
(e.g., Allen et al. 2009; Gnat & Sternberg 2009) and 
in turbulent mixing layers (e.g., Slavin et al. 1993). 

When departures from CIE are significant, the 
cooling efficiencies are no longer a function of just 
the gas temperature and metallicity, but instead de- 
pend on the specific time-dependent ion fractions. 
The use of convenient tables with known CIE cool- 
ing efficiencies must be replaced with a detailed 
computation of the non-equilibrium cooling, tak- 
ing into account all the relevant microphysical pro- 
cesses which include numerous emission lines, ther- 
mal bremsstrahlung, and ionization and recombina- 
tion cooling. This both requires the collection of 
a large set of atomic data for all the relevant pro- 
cesses, and is computationally complex compared 
with using look-up tables. 

Recently, first attempts have been made at in- 
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eluding nonequilibrium ionization physics in large 
scale hydrodynamical simulations, with applications 
for both eosmological metal-absorption through the 
warm/hot intergalactic medium (e.g., Cen & Os- 
triker 2006; Cen & Fang 2006) and emission- 
lines from galaxy clusters (e.g., Akahori & 
Yoshikawa 2010). Because of the complexity of 
nonequilibrium cooling, such simulation have so far 
not included self-consistent nonequilibrium cooling 
efficiencies. It is simpler to include nonequilib- 
rium cooling in hydrodynamical simulations by us- 
ing look-up tables for the nonequilibrium cooling 
rate (e.g., Oppenheimer & Dave 2009). 

In this paper, we present ion-by-ion cooling ef- 
ficiencies. We list the cooling efficiencies for each 
ion of the first 30 elements (H— Zn) individually. 
We present results for gas temperatures between 
10'' and 10^ K, assuming optically-thin, low-density 
conditions. The sum of ionic cooling efficiencies, 
weighted by the nonequilibrium ion-densities, then 
provides an efficient-to-compute and self-consistent 
nonequilibrium cooling efficiency. We also list the 
elemental cooling efficiencies assuming CIE condi- 
tions. These can be used to construct CIE cooling 
efficiencies for non-solar abundance ratios, or to es- 
timate the remaining cooling due to elements not 
included in any time-dependent computation. 

This paper is not a new calculation of the cool- 
ing functions (c.f. Schure et al. 2009). We do not 
present any new atomic data. Instead, we list the 
current cooling functions included in Cloudy (ver. 
08.00) in convenient online tables that are useful in 
any numerical computation in which the ion abun- 
dances are not in photoionization or coUisional equi- 
librium. It is the first time that the individual ionic 
cooling efficiencies are listed in an accessible format. 
As we describe below, the new framework that we 
present here will be periodically maintained and up- 
dated, as improved atomic data becomes available. 

The outline of this paper is as follows. In Sec- 
tion 2, we describe the computational method. In 
Section 3, we present the results for the ion- by-ion 
cooling efficiencies, and for the element-by-element 
CIE cooling efficiencies. We summarize in Section 4. 

2. METHOD 

We used Cloudy (ver. 08.00, Ferland et al. 1998) 
to compute the ion-by-ion cooling efficiencies of 
every ion of the first thirty elements (H— Zn). 
All the cooling processes considered by Cloudy 
are described in detail in Osterbrock & Ferland 
(2006), and include collisional excitations followed 
by line emission, recombinations with ions, col- 



lisional ionizations, and thermal bremsstrahlung^. 
The electron cooling efficiency includes the removal 
of electron kinetic energy via recombinations with 
ions, collisional ionizations, collisional excitations 
followed by prompt line emissions, and thermal 
bremsstrahlung. Cloudy does not include the ion- 
ization potential energies as part of the total inter- 
nal energy, but instead follows the loss and gain of 
the electron kinetic energy only. Therefore, in the 
definition of the cooling (see Osterbrock & Ferland 
2006; Gnat & Sternberg 2007) the ionization poten- 
tial energy that is released as recombination radia- 
tion does not appear. Only the kinetic energy of the 
recombining electrons contributes to the cooling ef- 
ficiency. On the other hand, kinetic energy removed 
via collisional ionization is included in the cooling. 
If ionization potential energy is considered as part of 
the total internal energy, then collisional ionization 
does not lead to a net energy loss, since the kinetic 
energy removed is merely stored as potential energy. 
Either way of accounting for the energy losses leads 
to the same net (i.e. cooling minus heating) cooling 
efficiency. 

For each ionization state i of each element E, 
we constructed a series of models for different gas 
temperature between 10"* and 10^ K. Each model 
includes only hydrogen and the element E. The 
abundances of all other elements are set to zero. 
We set an electron density rie = 1 cm~'^, regard- 
less of composition and ionization state. We define 
the ionization states of the element E so that the 
fractional abundance of the species Ei is 1, and the 
abundances of all other ions are 0. We further set 
the abundance of element E to be 10^^ larger than 
that of Hydrogen, so that nn = 10~^^ cm~'^ and 
n{Ei) — 1 cm~^. Effectively, each such model con- 
tains only the species Ei and free electrons at the 
specified electron temperature.^ 

For each element, we also compute the cooling 
efficiencies assuming CIE ion fractions. In these 
models, we set the abundance of the element E 
to be lO'^^ larger than that of Hydrogen (so that 

Cloudy is available at http:/ /www.nublado.orgT] 
The code documentation, including a full 3e^ 
scription of all cooling processes is available at 
http: / / viewvc.nublado.org/index.cgi / tags/release / cOS.OO / - 
docs/'.''root=cloudy (and will be described in Ferland et al. 
2011, in preparation). To get the full list of references to the 
atomic data in this version, see instructions in section 13.5 
of Cloudy's third volume of documentation (Hazy3_08.pdf). 

^ We verify that the contribution of Hydrogen to each 
model is negligible, by changing the abundance of the rel- 
evant species from 10^^ times rajj to 10^" times njj, and ver- 
ifying that the results remain unaltered. 
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again nn — 10^^'"' cm^'^ and n{E) — 1 cm~'^), 
and we force an electron density = 1 cm^ . 
We allow Cloudy to compute the CIE ion frac- 
tions x{Ei). For each element E, the ion fractions 
x{Ei) — n{Ei)/n{E), must at all times satisfy, 

Y^x{E,) = l, (1) 

i 

where n{Ei) is the density (cm^'^) of ions in ioniza- 
tion stage i of element E, n{E) = uhAe, nn is the 
total hydrogen density, and Ae is the abundance 
of element E relative to hydrogen. These models 
yield the total cooling efficiency due to the CIE ion 
distribution of element E. 

The tables presented in this paper provide easy 
access to different components of the total cooling 
function computed in version 08.00 of Cloudy. The 
cooling function included in Cloudy is constantly 
being updated as improved atomic data become 
available. We will update this table, keeping its 
current format, to provide ready access to these fu- 
ture calculations. This way, codes that can parse 
the current tables can be easily updated as better 
atomic data become available. 

3. lON-BY-ION COOLING EFFICIENCIES 

We have carried out computations of the cooling 
efficiencies Ag ion(2^) for each ion of the first 30 el- 
ements, H— Zn (with atomic numbers 1 — 30), as 
a function of temperature. The results are listed in 
tabular form i n Table [H The full table is avail able at 
http://wise-obs.tau.ac.i1/~orlyg/ion_by_ion/i and 
is divided into lettered parts A- AD, as is out- 
lined in Table [51 For each element E with atomic 
number Z , the first column in Table [T] lists the 
temperature, and the next Z + I columns list 
the cooling efficiencies, Ae.ion (erg s^^ cm"^) for 
the different ionization states, starting with the 
neutral atom, and ending with the fully stripped 
ion. The cooling rate per unit volume due to 
ionization state i of element i?, is then given by 
n(ion) TT-e Ae,ion (erg s~^ cm"'^). For example. Ta- 
ble UK shows that the cooling rate due to neutral 
hydrogen at a temperature of 1.1 x 10^ K is 1.33 x 
10~^^ n(H°) He erg s~^ cm~'^ whereas the cooling 
efficiency due to ionized hydrogen at the same tem- 
perature is 6.47 X 10^^^n(H+)ne erg s~^ cm~^. 

Figure [T] shows the ion- by- ion cooling efficiencies 
for hydrogen (upper panel) and for helium (lower 
panel). This Figure confirms that if neutral hy- 
drogen exists at high temperatures, it cools orders 
of magnitude more efficiently than ionized hydro- 
gen. For example, at a temperature of 10^ K, 



Table 1. Ion-by-Ion Cooling Efficiencies 



T 

(K) 


A,,Hi(H I) 
(erg cm'* s^-*-) 


Ae,Hn{H II) 
(erg cm'' 


A^,h(H at CIE) 
(erg cm"' s~^) 


1.00 X 10* 
1.05 X 10* 
1.10 X 10* 


4.59 X 10-24 
7.93 X 10-2* 
1.33 X 10-23 


6.26 X 10-25 
6.37 X 10-25 
6.47 X 10-25 


4.58 X 10-2* 
7.90 X 10-2* 
1.32 X 10-23 



Not e. — T able \T\ is available in i ts entirety at 
http://wise-obs.tau.ac.il/~orlyg/ion_byJon71 A portion is 
shown here tor guidance regarding its term and content. The 
full table lists the ion-by-ion cooling efficiencies for all the 
different ions of the first 30 elements (H— Zn), and for the 
elemental cooling efficiencies assuming CIE. For a guide see 
Table [2] 



Table 2. Cooling Data and Solar Elemental Abundances 



z 


Element 


Table 


Abundance (X/H)q 


Reference 1 


1 


Hydrogen 


lA 


1 




2 


Helium 


IB 


8.33 X 10-2 


B 


3 


Lithium 


IC 


2.04 X 10-8 


C 


4 


Beryllium 


ID 


2.63 X 10-" 


C 


5 


Boron 


IE 


6.17 X 10-10 


C 


6 


Carbon 


IF 


2.45 X 10-* 


A 


7 


Nitrogen 


IG 


6.03 X 10-5 


A 


8 


Oxygen 


IH 


4.57 X 10-* 


A 


9 


Fluorine 


11 


3.02 X 10-* 


C 


10 


Neon 


IJ 


1.95 X 10-* 


DT 


11 


Sodium 


IK 


2.14 X 10-"5 


C 


12 


Magnesium 


IL 


3.39 X 10-5 


A 


13 


Aluminum 


IM 


2.95 X 10-6 


C 


14 


Silicon 


IN 


3.24 X 10-5 


A 


15 


Phosphorus 


lO 


3.20 X 10-^ 


C 


16 


Sulfur 


IP 


1.38 X 10-5 


A 


17 


Chlorine 


IQ 


1.91 X 10-'' 


C 


18 


Argon 


IR 


2.51 X 10-'^ 


C 


19 


Potassium 


IS 


1.32 X lO-'' 


c 


20 


Calcium 


IT 


2.29 X 10-6 


c 


21 


Scandium 


lU 


1.48 X lO-'^ 


c 


22 


Titanium 


IV 


1.05 X lO-'' 


c 


23 


Vanadium 


IW 


1.08 X 10-* 


c 


24 


Chromium 


IX 


4.68 X lO-'' 


c 


25 


Manganese 


lY 


2.88 X 10-^ 


c 


26 


Iron 


IZ 


2.82 X 10-5 


A 


27 


Cobalt 


lAA 


8.32 X 10-* 


c 


28 


Nickel 


lAB 


1.78 X 10-6 


c 


29 


Copper 


lAC 


1.62 X 10-* 


c 


30 


Zinc 


IAD 


3.98 X 10-* 


c 



Note. — (1) References: A: Asplund et al (2005); B: 
Ballantyne et al. (2000); C: adopted from Cloudy, based on 
Grevesse & Sauval (1998); DT: Drake & Testa (2005). 
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neutral hydrogen has a coohng efficiency of ^ 
10""'^^ erg cm^ s~^, due to collisional ionizations and 
Lya emission. Ionized hydrogen has a cooling effi- 
ciency of ~ 2 X 10~^^ erg cm-^ s""'^, due to thermal 
bremsstrahlung emission. 




Fig. 1. — lon-by-ion Cooling efHciencies for hydrogen ions 
(upper panel) and for helium ions (lower panel). 

The total cooling due to a specific element de- 
pends of the ion abundances. Figure [2] shows the 
cooling efficiencies vs. temperature for the different 
carbon ions (see Table [T|^). The top panel shows the 
ion- by- ion cooling efficiencies, Ae^on- These can be 
used to compute the total carbon cooling efficiency 
for any composition. As an example, in the middle 
panel we show the carbon CIE ion fractions. The 
bottom panel then shows the contribution of each 
ion to the CIE carbon cooling efficiency Ae^c(ion) = 
a;ionAe,ion- For example, at a given carbon density 
n(C), the contribution of C'^"'" to the CIE cooling 



rate per volume is n(C) Ue Ae,c(C'^+) erg s~^ cm~'^. 
The sum of CIE cooling efficiencies over all carbon 
ions is shown by the thick gray curve in the lower 
panel. The carbon CIE cooling efficiency has two 
peaks. The first peak, at 6 x lO'* — 2 x 10^ K is due 
to C+, C^"*", and C'^"'". These ions are responsible 
for the familiar carbon peak in the solar-metallicity 
CIE cooling curve at 10^ K (see Section[3T|). The 
second peak, at ^ 10^ K, is due to C^+ and C^+, 
and is two order of magnitude lower. 

3.1. Element-by- Element CIE cooling 

The last column in each part (A- AD) of Table [T] 
lists the total CIE cooling efficiency of each element 
as a function of temperature. Figure [3] shows the 
CIE cooling efficiencies of the major coolants as a 
function of temperature, assuming a solar metal- 
licity. In making Figure [3l we adopt the elemen- 
tal abundances for C, N, O, Mg, Si, S, and Fe re- 
ported by Asplund et al. (2005) for the photosphere 
of the Sun, and the enhanced Ne abundance recom- 
mended by Drake & Testa (2005). For the other ele- 
ments we use the abundances reported by Grevesse 
& Sauval (1998). We list these abundances in Ta- 
ble p 

Figure |3] shows the familiar peaks in the CIE cool- 
ing efficiency due to different elements. The low- 
temperature peak at ~ 2 X 10^ K is mainly due 
to hydrogen Lya cooling. As the hydrogen neu- 
tral fraction becomes small, the contribution of hy- 
drogen Lya decreases. This peak is followed by 
peaks at 10^ 3 x 10^ 5 x 10^, and 1.5 x 10^ K 
due, respectively, to contributions of carbon, oxy- 
gen, neon, and iron. A second iron peak can be 
seen at ^ 10^ K. At higher temperatures cooling is 
dominated by thermal bremsstrahlung due to fully 
stripped ions. Contributions by other cooling el- 
ements are also shown in figure [31 For example, 
cooling due to helium peaks at a temperature of 
~ 8 x lO'' K. For solar metallicity gas the helium 
contribution to the total cooling is small compared 
with that of metal-line cooling. However, at subso- 
lar metallicities the relative contribution of helium 
is larger. Nitrogen, magnesium, silicon, and sulfur 
also contribute to the cooling below ~ 5 x 10^ K. 
The results for the elemental CIE cooling efficiencies 
are in qualitative agreement with previous compu- 
tations (e.g., Sutherland & Dopita 1993). Difference 
in detail are mainly due to differences in assumed 
atomic data (see Gnat & Sternberg 2007), affecting 
both the CIE ion fractions and the cooling efficien- 
cies of specific coohng processes. 

The upper thick gray curve shows the total con- 
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Fig. 2. — Cooling efficiencies of carbon ions. Upper panel: 
lon-by-ion cooling efficiencies for carbon ions. The cooling 
rate per unit volume is given by nen[an^e,ion{^on). Middle 
Panel: Carbon CIE ion fractions. Lower Panel: Carbon CIE 
cooling efficiencies for carbon ions, i^ion^e.ion- The cooling 
rate per unit volume is given by riencAe c(ion). The sum 
of CIE cooling efficiencies by all carbon ions is shown by the 
thick gray curve. 



tribution of all elements to the CIE, solar metal- 
licity, cooling efficiency. These results are identical 
to the cooling efficiencies computed by Cloudy (ver. 
08.00) assuming CIE and a solar composition (in- 
cluding all elements). For comparison, the dotted 
curve shows the CIE cooling efficiencies of Gnat & 
Sternberg (2007), which were computed using the 
cooling function included in Cloudy ver. 06.02. The 
agreement is excellent for T < 2 x 10® K, but some 
differences appear at higher temperatures. This is 
due to difference in the input atomic data between 
the two Cloudy versions. 

4. SUMMARY 

In this paper, we present computations of the 
cooling efficiencies of each ion of the first 30 el- 
ements (hydrogen-zinc) individually. We use the 
cooling functions included in Cloudy (ver. 08.00) 
to compute the cooling efficiencies as a function of 
temperature, between lO'* and 10® K, assuming op- 
tically thin conditions. 

The results are listed in tabular form in Table [1] 
(Section [3]), and are available in convenient online 
format through the electronic edition of the journal 
(for a guide see Table [2]) . For each ion, we list the 
cooling efficiency Aejon(r) (erg cm'^ s~^) as a func- 
tion of temperature. The total cooling rate for any 
ionic composition can then be computed by mul- 
tiplying the ionic efficiencies by the ion densities, 

'^e Sion '^(i0Il)^e,ion (erg S^^ Cm"^). 

As opposed to gas in CIE, for which the cooling 
efficiencies depend only on the gas temperature and 
mctallicity, for nonequilibrium conditions the cool- 
ing efficiencies must be evaluated locally depend- 
ing on the nonequilibrium ion fractions. A self- 
consistent computation therefore requires the col- 
lection of a large set of atomic data for all the rele- 
vant microphysical cooling processes, including nu- 
merous emission lines, thermal bremsstrahlung, and 
ionization and recombination processes. The results 
presented in this paper allows for an efficient esti- 
mate of the total cooling efficiency regardless of the 
ionization state. 

The tables presented here use the current atomic 
data set within Cloudy ver. 08.00. The tables pro- 
vide a flexible way to access the cooling of individ- 
ual species. The Cloudy atomic database is con- 
tinuously updated and new versions of these tables, 
using the same format, will be created as the atomic 
data are improved®. These tables will then provide 

^ Cloudy provides a means to generate its current atomic 
physics bibliography, as described in Section 2 above. 
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Fig. 3. — Element-by-element cooling efficiencies assuming CIE ion fractions and solar elemental abundance ratios (see 
Table [2]|. The total CIE cooling efficiency due to all elements is shown by the upper thick gray curve. The CIE cooling 
efficiency of Gnat & Sternberg (2007), which relied on Cloudy version 06.02 is shown by the upper dotted line for comparison. 
The differences between the two are due to updated atomic data included in Cloudy version 08.00. 
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easy access to future improvements as tliey occur. 

Tfie online tables are useful when constructing 
theoretical models in which nonequilibrium ioniza- 
tion plays a significant role, and can be used, for 
example, in models for radiatively cooling gas, con- 
duction fronts, fast shock waves, and turbulent mix- 
ing layers. They can also simplify the inclusion of 
self-consistent nonequilibrium cooling in large scale 
cosniological hydrodynamical simulations. 

In Section 13. 1[ we present the elemental cooling 
efficiencies as a function of temperature, for each of 
the first 30 elements (H— Zn) assuming CIE condi- 
tions. These results can be used to easily construct 
CIE cooling efficiencies for non-solar abundance ra- 
tios, as well as to estimate the cooling by elements 
not included in any time-dependent, nonequilibrium 



computations. 
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